SU(4) spin-orbit critical state in one dimension 
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Effect of quantum fluctuations concerned with the orbital degrees of freedom is discussed for 
the model with SU(4) symmetry in one dimension. An effective Hamiltonian is derived from the 
orbitally degenerate Hubbard model at quarter filling. This model is equivalent to the Bethe soluble 
SU(4) exchange model. Quantum numbers of the ground state and the lowest branch of excitations 
are determined. The spin-spin correlation functions are obtained numerically by the density matrix 
renormalization group method. It shows a power-law decay with oscillations of the period of four 
sites. The period originates from the interference between the spin and orbital degrees of freedom. 
The exponent of the power-law decay estimated from the finite size data is consistent with the 
prediction by the conformal field theory. 

Schultz-Mattis theorem, |5|,^| which states that the half- 
integer-S 1 spin chain, which has the translational and ro- 
tational symmetries, either has a singlet ground state 
with gapless excitations or has a finite gap with degener- 
ate ground states, corresponding to spontaneous break- 
ing of the parity. 

In the present paper we study an effective model of an 
orbitally degenerate Hubbard model in one dimension. 
By the density matrix renormalization group (DMRG) 
method and exact diagonalization (Lanczos method), we 
find a quantum critical state at the SU(4) symmetric 
point, which originates from the strong interplay between 
spin and orbital quantum fluctuations in one dimension. 



I. INTRODUCTION 

Recently the role of the orbital degrees of freedom in 
strongly correlated electron systems is attracting grow- 
ing interest. The increase of this attention is stimulated 
by the progress in the experimental studies of transition 
metal and rare earth compounds such as LaMn03, CeBg, 
and TmTe, which show various interesting properties as- 
sociated with the orbital degrees of freedom. 

In the 1970's Kugel and Kohmskii Q and Inagaki 
H studied an orbitally degenerate model to understand 
the magnetic structures of transition metal compounds 
within the mean field theory. They concluded that 
if orbitals ordered antiferromagnetically, then spins or- 
dered ferromagnetically, and vice versa. Recently Shiina, 
Shiba, and Thalmeier || have studied similar models in 
connection with a quadrupolar ordering of CeB@ and dis- 
cussed the phase diagram under an external magnetic 
field neglecting quantum fluctuations. 

In the case of LaMnC>3, the orbital ordering tempera- 
ture, To (~ 775 K), is much higher than the Neel tem- 
perature, Tn (~ 141 K), so the mean field theoretical 
approaches are considered to be a good starting point. 
On the other hand, for CeBg, To (~ 3.4 K) is the same 
order as TV (~ 2.3 K) and thus the interplay between 
spin and orbital quantum fluctuations may be important. 
Therefore it is necessary to consider the effects of quan- 
tum fluctuations more seriously beyond the mean field 
theory. 

Before considering the effects of the orbital degrees of 
freedom, we briefly summarize the properties of the one- 
dimensional single orbital Hubbard model for compari- 
son. In the limit of strong correlation at half filling, the 
model is reduced to the spin-1/2 antiferromagnetic (AF) 
Heisenberg model with SU(2) symmetry. This model is 
well known as a typical quantum critical system. The 
ground state of this model is singlet ([l 2 ] in Young's dia- 
gram representation) and the elementary excitations are 
gapless and triplet ([2 1 ]), so-called des Cloizeaux- Pearson 
modes [Q. These results are consistent with the Lieb- 



II. MODEL 

We start from the one-dimensional orbitally twofold 
degenerate Hubbard model with Hund rule coupling be- 
tween the two orbitals at the same site. This is the sim- 
plest model which possesses orbital degrees of freedom. 
Hamiltonian of this model is given by 

H = Ht + Hu + Hj 

= ( — f M+l c L<T C »+l a' a + H.cJ 

iota' a 

^ ^ ^J^iaa^ia' a' (1 0~aa' &aa' ) } 



iota aa 



(i) 



where c| (c, a(T ) denotes an electron creation (annihila- 
tion) operator with orbital a(=l,2) and spin a at the ith 



site, and 



is c , 



Si a denotes electron spin op- 



erator with orbital a at the ith site. Concerning the hop- 
ping matrix elements, the nearest neighbor hopping be- 
tween the same type of orbitals is assumed, tff +1 = t5 aa > . 
The simplest system which shows this property is illus- 
trated in Fig. [l]: the p x and p y orbitals along a chain par- 
allel to the z axis. We are interested in the case where t, 
U , and J are positive. 
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To study the region of strong correlation, we consider 
the limit of U, J 3> t at quarter filling. In this case charge 
degrees of freedom are suppressed and the system be- 
comes a Mott insulator. The effective Hamiltonian ob- 
tained by the usual second order perturbation is 




are the spin operators and 

Ti — — \Ci acr T aa > Ci a ' a J 

acta' 

are the pseudospin operators which describe the orbital 
degrees of freedom. In the above equations f are the 
Pauli spin matrixes. 

As a first step, we consider the case with the highest 
symmetry by taking the J — > limit. Then H c g becomes, 
neglecting a constant term, 

tt _ jcST p( S=1 /2) p(T=l/2) ,„s 

i 

where K = 2t 2 /U, P^ /2) = 2S % ■ S l+1 +l/2, and 

P^ /2) - ^ • f i+1+ l/2. P^ and P^ are 
the spin- 1/2 and the pseudospin- 1/2 exchange operators 
between the ith and (i+l)-th sites, respectively. 

Since the Hamiltonian (||) exchanges both S and T 
spins at the same time, the spin and orbital degrees of 
freedoms are combined into the SU(4) spin (denoted by 

~*( 3 / 2 ") 

Sp'') and H ee is described by using spin-3/2 exchange 
operators as follows: 

H sv{i) =Kj2H-:i 3/2 \ (4) 

i 

where 

p(S=3/2) _ 2/3(3/2) 3(3/2)^3 , 11/ 3( 3 /2) 3(3/2)^2 

_ 9/3(3/2) c(3/2)n_67 
8 Wt ' ) 32 • 

The Hamiltonian clearly has the SU(4) symmetry. We 
call this Hamiltonian the SU(4) exchange Hamiltonian. 

The exact ground state energy and the dispersion re- 
lations of the SU(4) exchange Hamiltonian have been al- 
ready obtained by the application of the Bethe ansatz 



technique to the higher spin-chain problems J7[. In this 
paper, we investigate this model as the coupled spin and 
orbital system with the strongest orbital quantum fluctu- 
ations. It is worth noting that the SU(4) exchange Hamil- 
tonian may play a similar role as the SU(2) AF Heisen- 
berg model for the single orbital Hubbard model. The 
assumptions that the hoppings of electrons are possible 
only between the same orbitals and the vanishing J pro- 
duce this SU(4) symmetry, independently of the strength 
of the Coulomb repulsion U . Generally speaking, in real 
materials the Hund rule coupling J is not small. How- 
ever, an understanding of the most symmetric case will 
be important for future studies of less symmetric cases 
corresponding to a finite J. 



. z-axis 




FIG. 1. A one-dimensional Hubbard model with twofold 
degenerate orbitals at each site. 



III. GROUND STATE AND EXCITATIONS 

To understand the physics of the present model, let 
us consider S and T spins as classical spins. The spin 
configurations where every adjacent two, either S or T, 
spins point the opposite direction have the lowest en- 
ergy. Thus the classical ground state energy is zero and 
the degeneracy of the ground states is macroscopic in the 
classical theory. This situation is different from usual or- 
bital and/or spin orderings discussed so far. Thus it is 
essential to examine the properties of Hamiltonian (Q) by 
unbiased methods. 

First, we calculate the ground state energy (E g . s .) 
by the DMRG method || in the subspace of 
(5 t 2 ot , T* t )=(0,0), where 5f oi = Sf and Tf ot = £\ T?. 
We take K — 1 as the energy unit and use here the 
open boundary conditions (OBC) to obtain sufficient ac- 
curacy by the DMRG method. The obtained results are 
shown in Fig.||, which shows that the ground state en- 
ergy per site (E g . s jN) and the surface energy are equal 
to —0.825 and 0.35(3), respectively. This ground state 
energy is, of course, consistent with that obtained by the 
Bethe ansatz, E g . s ./N ~ -0.825 12 §. 
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FIG. 2. Ground state energy as a function of N. The nu- 
merical errors, which are estimated from the truncation errors 
in the DMRG calculation, are less than 0.01%. The broken 
lines represent the linear fitting; E g . s . = —0.825 x iV + 0.35(3). 
The inset shows the excitation gap (A) as a function of 1/JV. 
Error bars are estimated from the truncation errors in the 
DMRG calculation. The dotted line represents the linear fit- 
ting, A = 7.3/N 



In Fig. H the ground state energies are plotted only for 
N = An, for which the minimum energy in the subspace 
of (Sf ot ,T* ot )=(l,l) is different from E g _ s ,. From Table I 
it is seen that the subspace (S^ ot ,T^ ot )=(0,0) is included 
in every irreducible representation, but (S* ot , T t * t )=(l,l) 
belongs to any irreducible representation except for [l 4 ]. 
Thus it is concluded that the ground state belongs to the 
[l 4 ] irreducible representation in the Young's diagram no- 
tation. Similarly, by calculating the ground state energies 
with changing (S*f ot , T* ot ), it is found that for N = An + 2 
the ground states belong to either [2 3 ] or [2 2 1 2 ], which 
are degenerate 10- and 6-fold, respectively (see Table I). 

These quantum numbers can be understood from the 
point of view of maximum antisymmetrization. That 
is, the irreducible representations thus obtained for the 
ground states are compatible with the simple fact that 
the more antisymmetric part one irreducible representa- 
tion has, the lower is its ground state energy in the sub- 
space, because the Hamiltonian (|]) is the sum of SU(4) 
exchange operators. To avoid complications coming from 
the degenerate ground states, we consider the systems of 
N = An in the following. In this case the ground state be- 
longs to the [l 4 ] and is a singlet. Since the Lieb-Schultz- 
Mattis theorem applies to this model, the excitations are 
expected to be gaplcss, provided that no other symmetry 
is broken, in the same way as the spin- 1/2 AF Heisenberg 
model with SU(2) symmetry. 

To estimate the excitation energy, we calculate the 
ground state energy and the first excited state energy by 
using the DMRG method for N > 20 and by the exact 
diagonalization (Lanczos method) for N < 16. We deter- 
mine the ground state energy and the first excited energy 
by the minimum energies of the states whose quantum 
numbers (S^ ot ,T^ ot ) are specified: (0,0), (1,1), etc. It is 



found that the first excited state belongs to [2 1 ! 2 ]. 

Though the DMRG is more suitable for OBC than 
periodic boundary conditions (PBC), here we apply the 
PBC in order to study the properties in the bulk limit. 
When we use the OBC, we get lower excitation energies 
than those shown in the inset of Fig.|2| but they corre- 
spond to excitations at the surfaces rather than those in 
bulk. From the inset of Fig. || we can conclude that the 
excitation gap (A) goes to zero as A ~ 7.3/N. 

In order to examine the properties of the excitations in 
more detail, we calculate the dispersion relation by using 
the Lanczos method with the use of translational sym- 
metry for the systems with the PBC. Figure || shows that 
the excitation spectrum has a "bactrian camel" structure 
and shows softening at q — tt/2. This structure is also 
known by the Bethe ansatz results . Corresponding to 
the softening at q = tt/2, the correlation functions would 
show a characteristic feature, namely, oscillatory behav- 
iors with a period of four, which we will discuss in the 
next section. 

Figure || shows that height of the left hump is always 
lower than that of the right one. To consider a possi- 
ble reason, we determine the irreducible representation 
of each state for N — 8 and 12. Quantum numbers as- 
signed for each q point in the dispersion curves are shown 
in Fig. H by the SU(4) Young's diagram representations. 
In these finite size calculations, the state at q = tt/2 and 
left part of the two humps always belongs to [2H 2 ] and 
the right one to [3 1 ! 1 ]. From Table I, the first excited 
states at q — tt/2 consist of the coupled spin and or- 
bital excitations in addition to the pure spin and orbital 
excitations and have the 15 x 2-fold degeneracy in total. 
The difference of the height of the two humps may be 
attributed to the difference of the irreducible representa- 
tions of the two humps. In fact, the left (lower) hump be- 
longs to the irreducible representation which has a more 
antisymmetric part than that of the right (higher) one. 
In the bulk limit, however, we expect that the two parts, 
< q < tt/2 and 7r/2 < q < tt, converge to the same dis- 
persion relation as is known by the Bethe ansatz solution 

TABLE I. Irreducible representations of the SU(4) sym- 
metry given by Young's diagrams (YD) and their relations to 
(S t ot,Ttot) and sfj t 2) representations [9], where S to t, T to t, and 
Stot are the magnitudes of the Stot = J2i = J2i ^* > 



'(3/2) 



and S. 



of each representations. 



(3/2) 



respectively. Here v is the degeneracy 



YD y (Stot, T tot ) Sfj t 2) 

[i 4 ] i W) 6 

[2 1 ! 2 ] 15 (0,1)9(1,0)9(1,1) 19203 

[2 2 ] 20 (0,0)9(1,1)0(0,2)9(2,0) 09 2 9 2 94 

[3 1 ! 1 ] 45 (0,1)9(1,0)9(1,1)9(1,2)9(2,1) 1919293939495 
[4 1 ] 35 (0,0)9(1,1)9(2,2) 092939496 



[2 2 L] 6 
[2 3 ] 10 



(1,0)9(0,1) 
(0,0)9(1,1) 



092 
193 
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FIG. 3. Dispersion relations. The symbols A, □, and O 
represent data for N = 8, 12, and 16, respectively. Young's 
diagrams (YD) show the irreducible representations. The 
numbers accompanied by YD show the degeneracy and x2 
represents the same weight from q and —q. The solid line 
represents the Bethe ansatz result [7]. 



IV. CORRELATION FUNCTIONS 

Now we move on to the behaviors of the correlation 
functions, (Sf ■ Sf + j) g . s ., where (• • -) g , s , denotes expecta- 
tion values for the ground state. Since Hamiltonian (||) 
has rotational symmetry with respect to both S and T 
spins, we consider only z components of spins. 

We use the OBC to get better accuracy in the DMRG 
calculations, but in this case we must keep in mind 
that the data contain the effects from boundaries. (Sf ■ 



S 



i+j/g-s. 



shows an oscillatory behavior with a period of 
four as a function of j. But the correlation functions also 
vary with a period of four with respect to i. That is, 
( s f ■ Sf +:j )g. a . is equal to (Sf +4 ■ Sf +4+j ) g . s . for any 

This behavior is caused by the standing wave with 
a period of four originating from the open boundaries. 
To remove such an artifact due to the OBC, we average 
(Sf ■ Sf + j)g. s , for one period with respect to i. Thus we 
define an approximate bulk correlation function as fol- 
lows: 



(Sf ■ Sf + j) 



'bulk 



1 3 

fc=0 



i+k ' ^i+k+j) g-s. 



(5) 



After this averaging procedure, we get the natural be- 
havior of the correlation functions as shown in Fig. |] . 

Because the results discussed in the previous section 
show that this model is gapless, we try to fit the enve- 
lope of (Sf • Sf + j)buik data with a power-law function ( 
j~ a ) by the least mean square method and get critical 
exponent a equal to 1.80 or 1.55 depending on using ei- 
ther the upper (j = 12, 16, and 20) or the lower (j = 14, 
18, and 22) data. We did not use the data of j = 24, 26 
and 28, because these sites are too close to the bound- 
ary. Due to finite size effects, the value of a depends on 



how to fit, but a is always between 1.5 and 2.0. From 
these results, we conclude that the asymptotic form of 
the correlation function is given by 



(Si ■ S i+ Abu. 



ilk 



COS (f j) 



a 



1.5-2.0, 



(6) 



in the bulk limit. 

Of course, (Sf ■ Sf + j)buik is equal to (Tf ■ Tf + j)b u iki 
because of the symmetry of Hamiltonian (^|) concern- 
ing the exchange between S and T. Furthermore, (Sf ■ 
Tf + j)buik always equals zero as is easily shown by the 
Wigner-Eckart theorem. In fact, the calculated values of 
(Sf ■ Tf + j)buik are almost zero, and the numerical errors 
for the values of (Sf ■ Sf + j)buik may be estimated from 
these values, which are less than 1% even at the farthest 
site from i. 




FIG. 4. Correlation functions ((Sf ■ Sf+j)bulk) for the sys- 
tems with N = 48 (+) and 60 (O). Broken and dotted lines 
are 0.144/j 1 ' 80 and -0.108/j 1 ' 55 , which are the results of the 
least mean square fitting using the upper and lower data, re- 
spectively. Numerical errors, which are estimated from the 
values of (Sf ■ T* + j)b u ik, are less than 1%. The inset shows 
the entire form of the correlation functions. The symbols x 
and □ represent data for N = 48 and 60, respectively, but 
they overlap nearly perfectly. 



Next we study the structure factor defined by 

JV/2-2 

S z (q)^ (Sf ■ Sf +3 )buik ■ . 

j=-N/2+3 



(7) 



As is seen in Fig. |[ S z (q) has a characteristic cusp struc- 
ture at q = 7r/2. This result is consistent with the soft- 
ening at q = 7r/2 in the dispersion relation. By Fourier 
transformation of Eq. @, the analytic form of S z (q) is 
given by 



S z (q)~S z (q=-) 



|? -if- 1 



2 r(a)sin- 



O q 



(8) 
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around q = ir/2, when 1 < a < 3. If a is greater than 
2, Eq. (|8|) does not show any cusp structure at q = tt/2. 
So a must be less than 2, since we clearly see the cusp 
structure of S z (q) which becomes sharper as the system 
size is increased. 

By the SU(4) conformal field theory, the critical expo- 
nent of the SU(4) spin correlation functions with q = it/ 2 
oscillations is obtained to be 3/2 flQ(| , which is consistent 
with the present numerical result. Although the corre- 
lation functions discussed in this paper are the .S'-spin 
correlation functions but not the SU(4) spin correlation 
functions, we can show that the exponent is the same for 
the two correlation functions. 




FIG. 5. Fourier transformation of the correlation function 
for the systems with N = 60 (□) and N = 48 (x). 

V. CONCLUSIONS AND DISCUSSIONS 

In conclusion, we have studied the quantum critical 
state for the coupled spin-orbit system. The quantum 
numbers of the ground state and the lowest branch of the 
excitations are determined. Furthermore, the spin-spin 
correlation functions are obtained explicitly for the first 
time by the DMRG method. It shows a power-law decay 
with a period of four, which originate from the interfer- 
ence between the spin and orbital degrees of freedom. 
The exponent of the asymptotic behavior is consistent 
with the prediction by the conformal field theory. 

In this paper we have investigated only the most sym- 
metric model, but it is more realistic to consider a model 
with lower symmetry corresponding to a finite J. In such 
a case the effective Hamiltonian is given by the S'-spin 
isotropic and T-spin Ising-type Hamiltonian (^), whose 
properties are not yet fully understood. 

Related to the SU(4) model, several models with lower 
symmetries have been studied Kawano and 

Takahashi jl2| discussed S'-spin isotropic and T-spin XY- 
type Hamiltonians to study the three-leg antiferromag- 
netic Heisenberg ladder and showed that such a model 



is gapfull and has exponentially decaying correlation 
functions. Kolczhuk and Mikeska [|l3) studied a spe- 
cial SU(2)xSU(2) symmetric Hamiltonian ^■(S i -S i+ i + 
3/4) (T • T + i + 3/4) and showed that this model is also 
gapfull. 

It may be possible to study the properties of these 
models in a unified way by introducing different type of 
anisotropies from the SU(4) symmetric point. For this 
purpose, it is highly desirable to develop an analytic the- 
ory around this symmetric point. 
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